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A SEMIEMPIRICAL MODEL FOR BOUNDARY LAYER DIFFUSION 
IN FORCED CONVECTION FLOW 
by Alex Vary 
Lewis Research Center 


SUMMARY 


This report pertains to mass transfer involving boundary layer convective diffusion 
as the rate-controlling factor. In this case, hydrodynamic and surface conditions and 
diffusing specie concentrations in the bulk stream influence the mass transfer rate 
through their effect on the concentration gradient in the boundary layer. Herein, semi- 
empirically based equations are used to evolve mass transfer relations that depend on 
the measurement of bulk stream concentration boundary conditions. 

It is observed that the instantaneous mass flux m is related to the axial concentra- 
tion gradient in the bulk liquid stream: nt = g'(dc/dx), (where c is the concentration 
and x is distance). The mass transfer coefficient g' is functionally related to the 
thermophysical properties of the fluid, hydrodynamic factors, and geometric properties 
of the reaction surface. Herein, g’ is expressed explicitly in terms of mean stream 
velocity and mass flow rate as well as in terms of dimensionless groupings of hydrody- 
namic and surface properties. 

The mass transfer coefficient is evaluated by means of the experimental conditions 
corresponding to the mathematical model m ^ = g f ( c L - Cq)/L, where m ^ is the inte- 

is the axial 


- L ^ ^ u 
grated mass flux (m/L)dx , L is the test section length, and 


C L c 0 


concentration change in an isothermal test section. Evaluation of g' then depends on 
computing mass flux m j__ per unit axial change of bulk stream concentration under con- 
ditions in which g' is axially invariant (or essentially invariant). 

Using the previously described method, g’ can be evaluated for a range of hydrody- 
namic and surface conditions to establish empirical values for the coefficients and expo- 
nents involved in functionally relating g T to hydrodynamic and surface parameters. The 
procedure requires actual measurements of concentration boundary conditions that have 
heretofore been circumvented. 



INTRODUCTION 


Corrosive attack and mass transfer effects are important factors to consider in the 
design and operation of energy conversion systems using liquid metals as working and 
heat exchange fluids. Although containment materials that are highly resistant to corro- 
sive attack by liquid metals have been identified and tested, the need to improve our 
understanding of corrosion and mass transfer processes remains. In high-temperature 
nuclear applications, more information concerning the influence of hydrodynamic factors 
on convective transfer of soluble materials is needed. Modes and rates of depletion and 
accumulation of alloy elements, interstitials, and reaction byproducts are of special 
interest in the development of reliable, long-life nuclear reactor systems cooled by liq- 
uid metals. 

Hydrodynamic factors will influence the rate of corrosive attack or deposit formation 
in liquid metal flows where boundary layer convective diffusion is rate controlling 
(refs. 1 and 2). The variation of convective mass flux with hydrodynamic factors de- 
pends on the relation between momentum and concentration distributions in the boundary 
layer. However, both theory and experimental experience are incomplete, and neither 
provides adequate information relating convective mass transfer to momentum transfer. 
Therefore, liquid metal investigators employ semiempirical relations based on analogies 
among mass, momentum, and heat transfer (refs. 3 to 6). But this approach requires 
the accumulation of much more experimental data than have thus far been obtained to 
verify these relations. 

Results of mass transfer experiments involving hydrodynamic variables in high- 
temperature liquid metals cannot always be unambiguously interpreted or extrapolated 
to predict long-term convective mass transfer results (ref. 7). This is partly due to 
complexities involved in conducting suitable experiments and measuring bulk stream con- 
centrations with high-temperature liquid metals under forced flow. Analytical and ex- 
perimental efforts have generally been directed toward circumventing actual concentra- 
tion measurements within test sections. The usual procedure also involves making 
semitheoretically based estimates of mass transfer coefficients. The resultant estimates 
of mass transfer coefficients are then used to predict mass flux or to compare results 
with theory. But this method depends on a priori assignment of values to numerical 
quantities appearing in certain empirical equations (e. g. , Sh = aRe n Sc m , the Sherwood 
equation). At best, predicted results can be expected to agree with measured mass 
fluxes only within an order of magnitude. 

This report is concerned with cases in which the aforementioned practice is unsatis- 
factory. Herein, the experimental techniques currently used with liquid metals are re- 
examined with a veiw toward securing more accurate convective diffusion mass transfer 
information. The major objective is to present convective mass transfer equations in a 
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form that will be useful in determining values for mass transfer coefficients empiri- 
cally. A mathematical model is proposed for the empirical evaluation of mass transfer 
coefficients and for facilitating the prediction of trends under various forced- flow bound- 
ary conditions. Application of the model would necessitate determining actual variations 
in axial concentration gradients. 

CONVECTIVE MASS TRANSFER EQUATIONS 
Preliminary Considerations 

This report deals with liquid metal corrosion in a conduit under forced flow wherein 
the rate -controlling process is convective diffusion in the liquid boundary layer. In this 
case, hydrodynamic variables and boundary conditions will influence the corrosion (or 
deposition) rates by their effect on convective mass flux from (or to) the conduit surface. 
The ensuing discussion will apply specifically to situations in which the diffusing sub- 
stance tends to saturate the liquid at the conduit surface or in which the diffusing sub- 
stance is precipitating from a supersaturated solution. The diffusing substance, as con- 
templated herein, may be the major constituent of the containment material, an alloy 
element, an interstitial specie, or a reaction byproduct. 

Hereinafter, for the purpose of the discussion, fully developed turbulent liquid flow 
through a circular conduit of uniform cross section is considered. The influence of 
boundary layer temperature gradients, heat flux, and concentration gradients of other 
components of a mixture on diffusion in the boundary layer are ignored. (It is assumed 
that the latter conditions do not necessarily limit the applicability of the equations given 
herein. ) 

The mass transfer boundary layer (i. e. , diffusion sublayer) is the domain within 
which the concentration gradient (the driving force) primarily resides. For convective 
mass transfer, both molecular and eddy diffusion affect the mass transfer rate in the dif- 
fusion sublayer. It is important to distinguish between the boundary layer associated 
with mass transfer and that associated with momentum transfer. Each can differ mark- 
edly from the other in both thickness and transfer profile. For turbulent flow, there will 
tend to be a large disparity between the momentum and concentration distributions near 
the mass transfer surface (ref. 8). 

Herein an equation is sought for the instantaneous steady-state mass flux at any 
given point or station axially along the conduit under fully developed, steady tempera- 
ture conditions. That is, temperature and, hence, thermophysical properties are con- 
stant at any given station. The local mass flux is presumed to depend only on hydrody- 
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namic variables (e. g. , velocity, mass flow, and fluid friction), concentrations, and the 
properties of the conduit surface. 

According to the Fick equation, the mass flux will be proportional to the local con- 
centration gradient in the diffusion sublayer: 


m - pD ( ^ \ (1) 

Wy=0 

where m is the diffusing mass flux, p is the bulk stream density, D is the diffusion 
coefficient, and 3c/9y is the concentration gradient normal to the reaction surface. 

(All symbols are defined in appendix A. ) 

The diffusion coefficient D consists of two terms under the turbulent conditions 
considered herein because eddy effects are superimposed on molecular diffusion in the 
sublayer (i. e. , D = D + D,). The molecular diffusion coefficient D is a function of 
concentration and temperature. The turbulent (eddy) diffusion coefficient Dj. is a func- 
tion of eddy velocity and characteristic eddy transport (mixing) length (ref. 8). There- 
fore, both D and 3c/3y will depend on the momentum transfer rate in the boundary 
layer, given that D^. is of the order of or greater than D m . 

Under the previously stated conditions, equation (1) is indeterminate unless the mass 
flux can be mathematically related to the momentum flux. However, both theory and ex- 
perimental experience are incomplete, and neither provides adequate information relat- 
ing convective mass transfer to momentum transfer (ref. 9). Therefore, one might rely 
on some assumed analogy between mass and momentum or energy transfer. Since rele- 
vant fluid friction and heat-transfer data have generally been collected for zero mass 
transfer, the validity of the aforementioned analogies will be taken as essentially limited 
to dilute concentrations of the diffusing substance. The amount of solute introduced is 
assumed to have a negligible effect on the density of the bulk liquid even at saturation. 
These conditions are considered to be satisfied in nuclear liquid metal systems in which 
the diffusing specie concentrations are generally of the order of 1000 ppm or less. 


Correlation of Mass Flux With Hydrodynamic Factors 

In this section, an explicit relation is developed for convective mass transfer in 
terms of the experimental operating variables, such as velocity and mass flow rate, and 
in terms of nondimensional parametric groupings. The resulting equations will serve to 
define a mass transfer coefficient in terms of the axial concentration gradient along the 
flow conduit. 

Consider first the widely accepted mass transfer relation based on the following em- 
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pirical heat-transfer relation for turbulent flow in round tubes (ref. 10): 

Nu = aRe n Pr m (2) 

The mass transfer analog of equation (2) is formed by substituting appropriate nondimen- 
sional mass transfer groupings for the Nusselt Nu and Prandtl Pr numbers: 

Sh - aRe n Sc m (3) 

The transformation of Pr to Sc is accomplished by substituting pD for k/Cp. The 
transformation of Nu to Sh is accomplished by substituting g/pD for h/k, where the 
mass transfer coefficient g is defined as follows (ref. 11): 

m = g(c w - c g ) (4) 


Hence, equation (3) becomes 


m - apDRe n Sc m 



(5) 


Since (c - c )/d is clearly an approximation for 3c/3y, the term (aJte n Sc m ) must rep- 
resent an adjustment factor, since m in equation (5) is necessarily identical to m in 
equation (1). Thus, equation (5) is simply a modified version of the Fick equation. 

It is advantageous to relate m to the bulk stream axial concentration gradient in- 
stead of to the radial gradient across the boundary layer. An equation that is analogous 
to equation (5) can be obtained to relate mass flux m to the axial concentration gradient 
dc/dx in the bulk stream; for example, 


m = 


M dc 
i\d dx 


(6) 


Under the simplifying assumption that the mass m ird dx is instantaneously distributed 

o 

in the differential volume (1/4)7 rd dx, the coefficient a might be taken equal to unity. 
According to the analysis of appendix B, however, a is a function of Re and Sc; hence, 

m = a. — Re n Sc m — (7) 

77 d dx 

(Note that the coefficient a and exponents m and n in equations (2), (3), and (7) do 
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not necessarily have the same numerical values. ) Mathematically eliminating d be- 
tween M fd and Re 11 in equation (7) results in the desired equation for m , which is 
explicit in velocity, mass flow rate and axial concentration gradient (note that pU = M/A 

o S 


and A = (1/4)7 Td ): 


m = a’Se m p (1 - n > /2 ^ n (MU)<“- 1 >/ 2 *= (8) 

dx 

The coefficient a' differs from a in equation (7) in that it is a function of the exponent 

th 

n, since a' includes numerical quantities raised to the n LI power. 

Equation (8) reduces to an equation for mass flux in terms of a mass transfer coef- 
ficient and the axial concentration gradient (note that the axial concentration change is 
not a driving force in the same sense as is the radial concentration change across the 
boundary layer): 


** = g’ ^ (9) 

dx 

This new mass transfer coefficient g' is expressed explicitly in terms of velocity and 
mass flow rate: 

g* = a’Sc m p< 1 ' n)/ V n (MU)< n+1 >/' 2 (10) 

More generally, the mass transfer coefficient can be expressed in terms of nondimen- 
sional groupings that include the Reynolds number (see appendix C): 

g r = bpi/Sc m Re n P e (11) 

Equation (11) is derived directly from equation (10), and the equations are mutually ex- 
clusive with respect to MU and Re. The dimensionless parameter P is related to the 
geometric properties of the corroding surface. As explained in appendix C, the param- 
eter P can be expressed in different ways (e.g. , e/«/, L /d, d^/d 2 > etc.). 

In either equation (10) or (11), g’ is seen to depend on thermophysical properties of 
the fluid and on hydrodynamic and surface geometry factors. Since surface factors are 
accounted for through P in equation (11), surface factors must be implicit in a’ in 
equation (10). Obviously, this means that equation (10) applies only to cases with iden- 
tical surface conditions, if a* is to be taken as constant. 
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APPLICATION AND INTERPRETATION OF ANALYTICAL RESULTS 


Equation (9) is the defining equation for the mass transfer coefficient g'. Equa- 
tions (10) and (11) are expressions for g T in terms of experimental variables and pa- 
rameters that influence its value under various conditions. In this section, methods for 
evaluating g’ in a forced-flow test section are considered. A mathematical model 
suitable for this purpose is proposed and discussed. Subsequently, the variation of g’ 
and, hence, m with hydrodynamic and surface factors is discussed. Particular atten- 
tion is given to the ambiguities in past liquid metal corrosion experiments. These ambigu- 
ities apparently resulted from circumvention of the concentration measurements which are 
required by the present model. 


Evaluation of Mass Transfer Coefficient 

In any mass transfer experiment involving a forced-flow circuit, one or more con- 
duit segments or test sections are usually considered. The test section should contain 
a suitable reaction surface with respect to which mass removal or deposition rates may 
be determined. It should also provide for meaningful measurement of concentration dif- 
ferences for some predetermined length. It is by these means that the empirical evalu- 
ations of the mass transfer coefficient being sought can be obtained. 

Therefore, equation (9) is adapted to a conduit segment of length L along which the 
hydrodynamic and surface factors are constant. Also, for simplicity, let it be given that 
the conduit segment is isothermal and, hence, that thermophysical properties are con- 
stant. Under these conditions, the mass transfer coefficient can be considered a con- 
stant; therefore, equation (9) will yield the following expression for the integrated mass 
flux along the conduit segment L: 

= J' ~ dx = ^ (C L " V (12) 

where Cq is the bulk concentration at x = 0 and c^ is the bulk concentration at x = L. 
Equation (12) defines g' in terms of quantities that can serve as a basis for correlating 
mass flux with the hydrodynamic parameters included in equation (10) or (11). 

Caution must be exercised in using equation (12), since the net measured mass flux 
is a time-averaged as well as a space-averaged quantity. As corrosive attack progress- 
es, the mass flux may change because of changes in the nature of the corroding surface. 
The rate-controlling process may cease being boundary layer convective diffusion because 
of chemical changes in the substrate or the liquid. Or, surface variations may produce 
hydrodynamic changes that alter the mass transfer coefficient (see appendix D). 
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When there are interim changes in the mass transfer rate, there will generally be 
corresponding changes in c^ - Cq. It is apparent from this and the previous considera- 
tions that the validity of mass transfer correlations will depend on ascertaining whether 
there is a variation of the axial concentration gradient with respect to time as well as 
with respect to hydrodynamic and surface conditions. By analogy, one should remember 
to consider concentration measurements to be as important in mass transfer as temper- 
ature measurements are in heat-transfer investigations. 

The previous discussion indicates that the mass transfer coefficient can be accurately 
evaluated only under rather special conditions, even if concentration measurements are 
made. Furthermore, it should be apparent that a proper evaluation of g' depends on 
computing mass flux per unit axial change of bulk stream concentration under conditions 
in which g f can be taken as being axially invariant. Hereinafter, some data are 
presented to indicate appropriate values for the exponent n in equations (8) and (10). 


Variation of Mass Flux With Reynolds Number 

One purpose of equations (9) to (12) is to provide a format for experimentally evalu- 
ating the exponents m, n, and e. However, it has been suggested that an appropriate 
value for n can be assigned on the basis of the heat-transfer analogy (ref. 1). The tur- 
bulent heat-transfer exponent of 0. 8 has been adopted by many investigators for corre- 
lating convective mass transfer with Reynolds number. Although this particular value 
is not unequivocally confirmed by experimental results obtained heretofore, we will 
tentatively adopt it herein. The uncertainty with regard to the Reynolds number expo- 
nent is illustrated by the case cited below. 

The authors of reference 12 obtained mass transfer data for stainless -steel speci- 
mens exposed to flowing lithium in a forced-circulation loop. There was mass removal 
from specimens in a hot-sample section and mass deposition in a cold-sample section. 
Data were obtained for Reynolds numbers between 1000 and 16 000. Both mass removal 
and deposition were patently influenced by hydrodynamic factors, which indicated that 
boundary layer convective diffusion was rate controlling in each instance. The deposi- 
tion rate increased linearly with Re®' ® for each of several temperatures from 328° to 
490° C in the cold-sample section. However, mass removal in the hot-sample section 
increased parabolically with Re®' ® for each of several temperatures from 510° to 
612° C. The mass removal data indicated a higher power of Re (e.g. , approximately 
twice 0. 8). Or, if 0. 8 is valid for the Reynolds number exponent, then other variables 
must have influenced the mass transfer rate. But this question remains open because 
reference 12 provides no data relative to concentration gradients in the test sections. 
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Moreover, there is insufficient information concerning the method used for varying Re, 
which was defined in terms of sample length (i. e. , it is not clear whether sample length, 
velocity, or both were variables). 


Variation of Mass Flux With Velocity 

Many investigators prefer to correlate mass flux directly with velocity instead of 
Reynolds number. A common practice has been to adopt the Reynolds number exponent 
n for the velocity. According to equation (8), however, the velocity exponent equals 
(n+l)/2. If n has values near 1, the difference between n and (n+l)/2 is small and 
either value will suffice for practical purposes. 

If mass flow rate M is held constant, equation (8) states that mass flux will vary in 
direct proportion with provided that dc/dx and other factors are constant. 

But note that, in the same conduit (or in conduits with identical diameters or cross sec- 
tions), the mass flow rate will be directly proportional to velocity (since M = pUA s , 

A = constant). In this case, equation (8) states that m varies directly with u( n+ *), 

provided that dc/dx and other factors are constant. If n is 1, then in the former case 

o 

m varies directly with U, and in the latter case m varies directly with U . 

Both approximately linear and parabolic variations of mass flux with respect to ve- 
locity have been observed. In one investigation involving sodium in forced-flow stainless- 
steel loops (ref. 13), the mass flow rate was held constant so that the only hydrodynamic 
variable was velocity. Data were obtained with a number of loops under the aforemen- 
tioned conditions, and an analysis showed that the mass removal rate varied approxi- 
mately linearly with velocity. The velocity exponent was computed to be 0. 884. This 
agrees quite well with the velocity exponent of 0. 9 obtained in equation (8) if n is as- 
signed the value of 0. 8 as is usually done when using equation (5). On the other hand, 
data reported in reference 14 indicate a parabolic variation of mass removal rate with 
velocity for nickel alloy conduits. According to equation (8), this result should be ex- 
pected if n is approximately 1 and M varies directly with U while dc/dx is fixed. 

It appears that the foregoing results are consistent with the assumption that varia- 
tions in the overall axial gradients were negligible relative to the variation in velocity. 
Otherwise, the results would not agree so well with those predicted by equation (8). But 

this would imply that changes in velocity produce compensating changes in c - o, 

w o 

This would be necessary in order that dc/dx remain essentially fixed (see eq. (B6)). 
However, neither reference 13 nor 14 provides the concentration information required 
to support assumptions concerning either the constancy or mode of variation of axial 
gradients. 
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Variation of Mass Flux With Mass Flow Rate 


Tests have also been performed to determine the effect of variations in mass flux 
with mass flow rate at some fixed velocity, but the results are inconclusive. Two tests 
specifically devised for this purpose are reported in references 15 and 16. Special cor- 
rosion cartridges were installed in forced-flow stainless-steel loops with sodium. The 
cartridges contained annular channels that produced different mass flow rates at a fixed 
velocity. The cartridge design in the loop of reference 15 differed from that in the loop 
of reference 16. In the former loop, the mass removal rate increased in direct propor- 
tion with mass flow rate. In the latter loop the mass removal rate was essentially unaf- 
fected by mass flow rates differing by a factor of 2. 

Equation (8) states that m should vary directly with M if U and dc/dx are con- 
stant, but variations of M might produce inverse variations of dc/dx. From equa- 
tion (8) it is then possible that under some conditions m will appear to be invariant with 
respect to M. In the two tests cited, the difference in cartridge geometries could prob- 
ably partly account for the disagreement in results. Unfortunately, not enough informa- 
tion exists to explain the manner in which dc/dx differed between the two cases. 


Variation of Mass Flux With Surface Factors 

Convective mass transfer involving eddy diffusion will be influenced by surface 
roughness and irregularities and the geometric character of the test section or conduit. 
These factors can influence mass flux independently of the nominal Reynolds number and 
other parameters included in the mass transfer coefficient g’. However, the exponen- 
tial relations of mass flux with respect to surface factors have been studied even less 
than has mass flux with respect to hydrodynamic factors. 

Surface geometry effects are represented by the P parameter introduced in equa- 
tion (11) (see appendix C). Even when not represented explicitly (e. g. , eq. (10)), sur- 
face geometry effects should nevertheless be considered factors in the mass transfer 
coefficient g’. The definition of P will vary because it depends on the surface and 
geometric parameters that characterize the test section type being investigated. More- 
over, the P parameter is subject to change if the surface conditions change with time 
because of mass removal or deposition (e. g. , enroughening). 

There are many instances where anticipated correlations with Re or Sc have not 
been realized because surface and geometry variables altered the expected convective 
mass transfer rate. The authors of reference 3, for example, found appreciably in- 
creased convective mass transfer rates due to enroughening of benzoic acid surfaces by 
the formation of small fissures. In liquid-metal corrosion loops, locations with cracks 
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and crevices show higher rates of attack (ref. 15). Also, higher rates of mass removal 
or deposition are found on leading (front) surfaces of specimens and in locations with 
flow disturbances (ref. 17). 

The parameter P might correspond to the so-called downstream effect or x/d 
term included in some definitions of the mass transfer coefficient. In references 4 
and 17, for example, this x/d ratio is simply a normalized length against which vari- 
ations in m, are referred. Such a term is unnecessary in an expression for g* if a 
mathematical model of the form of equation (12) is adopted. Therein, the axial variation 
of mass flux is already taken into account by the definition of as the integrated 

average mass flux. 


MEASUREMENT AND ANALYSIS OF CONCENTRATION VARIATIONS 

Various techniques have been employed to detect impurities and solute concentra- 
tions that influence corrosive mass transfer in liquid metals. These include methods for 
detection and analysis of dissolved or entrained species in the bulk liquid at ppm levels. 
Devices used include plugging indicators, conductivity meters, cold traps, cover gas 
chromatographs, and drawn-sample vacuum- distillation apparatuses, and so forth (ref. 17). 

However, techniques for in-line monitoring of concentration variations along test 
sections have not been sufficiently developed. There are two basic methods that could 
serve this purpose: (1) spaced (periodic) direct sampling and analysis of the bulk stream 
through openings into the test section conduit, and (2) continuous indirect sensing of 
relevant physical or radiological changes that correlate well with solute concentra- 
tions. New and novel experiment designs and advanced instrumentation would be re- 
quired. Previous results indicate that the needed probe, instrument, and analytical sen- 
sitivities can probably be realized. Suitable in-line techniques and devices are described 
and discussed in reference 7 for the ppm range. However, future needs will undoubtedly 
involve sensitivities in the ppb range. The prospects of suitable devices in this range 
are unknown. 

It has been suggested that concentration variations might be deduced from post-test 
inspection and analysis of mass transfer results. For example, the variations of pene- 
tration depth or deposit height with axial distance downstream of the test section inlet 
might be interpreted as an indication of similar variations in the bulk liquid concentra- 
tion during testing. However, this technique may be satisfactory only in special and 
rather simple situations. The inclusion of test specimens or geometric complexities 
may invalidate assumptions of simple correspondence between dc/dx and mass removal 
or deposition profiles. 

In practical applications and in experimental forced-flow circuits, there will be two 
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extreme temperatures producing thermal gradient mass transfer (refs. 1 and 2). Asso- 
ciated with each temperature is a known (or assumed) equilibrium concentration of the 
diffusing specie(s). This knowledge is often used to define thermally induced concentra- 
tion differences that give rise to mass transfer; hence, it is used to circumvent actual 
measurement of c L - Cq if one is attempting to evaluate the limiting case. For ex- 
ample, in mass transfer experiments, the difference between the equilibrium concentra- 
tion for each of the two extreme temperatures can be substituted for c^ - Cq to produce 
correlations between g' and m provided that the following assumptions are reasonable 
approximations for the real system. First, liquid entering the test section has a concen- 
tration corresponding to saturation at the lower temperature (higher temperature for 
deposition). Second, the test section is practically isothermal. Third, liquid leaving the 
test section has a concentration corresponding to saturation at the higher temperature 
(lower temperature for deposition). 

In any test situation meeting the aforementioned conditions, it is desirable to attain 
essentially equilibrium saturation in the liquid exactly at the end of the test section for 
all flow conditions. However, no a priori criterion is available for selecting suitable 
test section lengths for each flow condition. In lieu of making all the appropriate meas- 
urements, therefore, there is no totally accurate means of evaluating the mass transfer 
coefficient. 


CONCLUDING REMARKS 

Heretofore, liquid- metal corrosion experiments and analyses have circumvented 
making actual measurements of mass transfer coefficients and concentration boundary 
conditions. The mathematical model proposed herein requires that convective mass flux 
be computed with respect to the axial concentration change along the test section. It is 
suggested that otherwise it is futile to expect valid and consistent correlations between 
convective mass flux and hydrodynamic (velocity, mass flow rate, and Reynolds number) 
and surface (roughness and form drag) factors. This situation arises because variations 
in hydrodynamic and surface factors produce variations in both boundary layer and bulk 
stream concentration gradients. 

Obtaining the needed concentration measurements is difficult with high-temperature 
liquid-metal systems. However, it is worth considering the adoption of experimental 
techniques and devising instrumentation that could aid in acquiring the axial concentration 
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data required by the mathematical model proposed herein. If done, this should facilitate 
more precise prediction of mass transfer rates with respect to variations in hydrodynam- 
ic and surface factors. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 4, 1970, 

120-27. 
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APPENDIX A 


SYMBOLS 


[The symbols are defined in terms of the basic dimensions of mass, M; length, L; 
time, T; force, F; and temperature, 9. ] 


A 

s 

cross-sectional area of con- 


gravitational constant, ML/FT' 


duit, IT 

h 

heat-transfer coefficient, 




F/LT0 

a, a*, a M 

dimensionless coefficient 





K 

numerical coefficient, M/FT 

b 

dimensionless coefficient 





k 

thermal conductivity, F/T0 

c 

point (or bulk stream) concen- 





L 

test section axial length, L 


tration, M/M 



c p 

specific heat at constant pres- 

l 

specimen length, L 


sure, FL/MT 

M 

mass flow rate, M/T 

c w 

concentration at reaction sur- 

m 

numerical exponent 


face (wall), M/M 

m 

mass flux, M/L 2 T 

c s 

diffusion sublayer concentra- 

Nu 

Nusselt number, hd/k 


tion at y = 6 (or y > 6 




for eq. (5)), M/M 

n 

numerical exponent 

D 

2 

diffusion coefficient, L /T 

P 

dimensionless parameter 



Pr 

Prandtl number, c nj//k 

D m 

molecular diffusion coeffi- 




cient, L 2 /T 

Re 

Reynolds number, d U/v 



Sc 

Schmidt number, v/D 

D t 

turbulent diffusion coeffi- 



l 

cient, L 2 /T 

Sh 

Sherwood number, g d/pD 

d 

conduit diameter, L 

t 

time, T 

e 

numerical exponent 

U 

mean stream velocity, L/T 

f 

friction factor 

u 

point axial velocity, L/T 

g 

mass transfer coefficient, 

w 

point axial mass flow rate, 


M/L 2 T 


M/T 

g’ 

mass transfer coefficient, 




M/LT 




14 



x axial distance along conduit, L 

y radial (normal) distance from sur- 

face, L 

6 diffusion sublayer thickness, L 

e root -mean- square roughness 

height, L 


o 

v kinematic viscosity, L /T 

p 

v molecular viscosity, L /T 

2 

turbulent viscosity, L /T 

O 

p density of liquid, M/L 
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APPENDIX B 


RELATION OF MASS FLUX TO AXIAL CONCENTRATION GRADIENT 

The relation between the mass flux from the conduit wall and mean axial concentra- 
tion gradient in the liquid stream is derived in this appendix. 

Let us consider a segment of a circular conduit with diameter d as illustrated in 
sketch (a). This analysis will be concerned with the differential annular volume indi- 


m 


^Circular conduit ; 

r 

d 6 


Bulk flow 
y -IV 

( ) 

1 C 6 


c 6 +dc 6i . 

\y 


7%% 


| S. 


^-dx 


la) 


cated. The concentration at the reaction surface c w is constant for the entire conduit 

segment. The diffusion sublayer thickness is 6 and the concentration at y = 6 changes 

from c. to c. + dc. over dx. 

0 0 0 

The mass balance for a differential annular element within the diffusion sublayer is 


c dw+ 7 rdm dx = 


[c + — dxldw + 7T d 
\ dx 


m 


dm 

9y 


dy dx 


(Bl) 


where m is the mass flux through the diffusion sublayer into the bulk stream and dw 
is the mass flow rate in the differential annular element of the diffusion sublayer. The 
following expression for m is obtained from equation (Bl): 


= — f 6 ^ — dy (B2) 

Trd Jq dy 9x 


The variation of mass flow rate w with y in 5, is 


— = pu7r( d - y) ~ pvnrd, 0 < y < 5 

dy 


(B3) 
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where 6 « d. From reference 10, 

(*i\ =*s! (a 

Wy=0 2 " 

rY 2 , 

Let du/dy in the diffusion sublayer equal (du/dy)y_Q, integrate u = ^ (fU /2u)dy, 
and substitute for u in equation (B3). 

It will be assumed that dc/dx increases linearly with y in the diffusion sublayer, 
that is, 


dc = y ^6 (B5) 

3x 5 dx 

It will be further assumed that dc^/dx is related to the mean bulk concentration gradient 
dc/dx through the Reynolds and Schmidt numbers, that is, 

dc K Hr 

— ^ - <p(Re, Sc)— (B6) 

dx dx 

Combining equations (B2) to (B6) and integrating yield 

m = fi U - 2 j 5 - 2 cp{ Re, Sc)— (B7) 

Qv dx 

Since pU 3 = (4/7 t</ 3 )MU, we have 

m = 2/L <p(R e , Sc)— (B8) 

7 id 3 vd dx 


or, 


m = -M 2f /_6 V Re <p(R e , Sc)— (B9) 

77 3 \d / dx 

From reference 8, the relation between the diffusion sublayer thickness Q and the vis- 
cous sublayer thickness 6 q is 

6 = 6 0 Sc" 1/3 (BIO) 
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where 5 q > 5. If the viscous sublayer is assumed to be laminar, then according to ref- 
erence 10, 


-j - (5 Re -1 ) \ 
d f f 


(Bll) 


Combining equations (B9) to (Bll) yields 

m -M.100 Re _1 Sc“ 2/,3 ^(Re, Sc)^. (B12) 

ird 3 dx 


It is deemed reasonable to assume that cp(Re, Sc) is a simple power function of Re and 
Sc so that equation (B12) becomes 


„ M „ n c m dc 
m - a — Re Sc — 


(B13) 


(The exponent n probably ranges from 0 to 1, but we refrain from assigning any particu- 
lar value for n, or for a or m, at the present stage. ) 
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APPENDIX C 


SURFACE GEOMETRY PARAMETER 


Examination of equation (B9) shows that the friction factor f and, hence, surface 
roughness or form drag will influence convective mass transfer. Moreover, for turbu- 
lent flow over rough surfaces we should not expect equation (Bll) to hold since the vis- 
cous sublayer may not be laminar. In general, therefore, f will not cancel out as in 
equation (B12) and some power of f should be implicit in the coefficient a’ of equation (8). 

Consider the following equation which is equivalent to equation (8) except for f : 


m 


i”f e Sc m n( 1 - n)/2 ^V (n+1) /^Y n+1)/2 — 


dx 


(Cl) 


The coefficient a" differs from a' by the factor (ir/4)^ n+ ^^ and by the purely 
numerical coefficient assumed to be associated with f e . Let us proceed by noting that 
equation (Cl) is unaltered if combined with the factor (pU^ /pTJy n+ ^ ^ : 


m 


- a"f e sc”pY^Y n+1) ( ^ro_y ,l+1)/2 * 

\ V } \ A _ tt 2 / dx 


(C2) 




Since t/U/v equals Re and MU/A pU z equals unity, equation (C2) may be written as 


m = a"f e pi^Sc ra Re n+1 — <C3) 

dx 

The friction factor f may be a simple function of the Reynolds number as in the 
Blasius equation, f = 0.079/Re which applies to smooth circular tubes. On the 
other hand, f may be totally independent of Re as in the Hopf equation, f = 0. 01(2e/</)®'^^ 
which applies to turbulent flow in rough tubes (ref. 10). Clearly, the exponents of f and 
Re will change as surface conditions are altered. Furthermore, in the case of conduits 
with in-stream corrosion specimens or other objects, the Reynolds number is usually 
written in terms of the object's characteristic dimension, say l, so that Re = l\J/u. But 
we can retain Re in its conventional form by combining it with the factor l/d. 

In accordance with the preceding considerations, equation (C3) is rewritten in terms 
of the parameters Re, Sc, and P: 


I 
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(C4) 


m = bpi/Sc m Re n P e — 
dx 

The exponent n in equation (C4) is not necessarily the same as n + 1 in equation (C3). 
Likewise, the coefficient b in equation (C3) will differ from a" or a’ in all the pre- 
vious equations ((8), (10), (Cl), and (C2)). In any case, the exponents and coefficients 
should be experimentally determined. (See appendix D for further comments concerning 
the exponents. ) 

The parameter P is related to surface and conduit geometry particulars that char- 
acterize the test section or conduit segment of interest. The surface geometry factors 
involved in P are those which would affect the boundary layer concentration gradient 
and, hence, dc/dx. Therefore, P should be differently defined for each type of reaction 
surface according to roughness, presence of in-stream corrosion specimens, annular 
conduits, protrusions, corners, and so forth. We can thus account for mass flux varia- 
tions due to a variety of typical surface geometry factors encountered in practical and 
experimental situations. For example, given uniform conduits that differ only in surface 
roughness, P might be taken simply as the roughness factor e/d. If specimen length 
effects need to be taken into account, then P might be set equal to l/d. If annular 
channels are used, then diameter ratios d-^/d^ would be appropriate. Obviously, P 
may be composed of combinations of suitable terms. 
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APPENDIX D 


VARIATIONS IN MASS TRANSFER EXPONENTS 

After any corrosive process has been in progress for a sufficiently long period, the 
original surface may have been altered by particulate mass removal, leaching, cracking, 
deposition, and other like processes. Therefore, as the attack goes on, the mass trans- 
fer coefficient will change even though the rate- controlling process, velocity, mass flow 
rate, and thermophysical properties are unchanged. 

One mode by which the mass transfer coefficient varies with respect to surface con- 
ditions is due to the variabilities of the exponents m, n, and e in equation (Dl) with 
respect to surface conditions: 


g' = bpj/Sc m Re n P e (Dl) 

This situation is analogous to that found in heat-transfer studies where the roughening of 
an initially smooth surface alters the relation among Nu, Pr, and Re. This is reflected 
by significant changes in the empirical exponents m and n. Like changes should be ex- 
pected in the exponents involved in the convective mass transfer coefficient. For ex- 
ample, it can be shown that the Reynolds number exponent may tend to increase (e. g. , 
from 0. 8 to 1.0) as the conduit surface is enroughened by mass removal. For the ve- 
locity and mass flow rate product MU common to most liquid metal applications, even 
small changes (e. g. , 0. 1) in (n+l)/2 in equation (D2) can produce significant changes 
in the mass transfer coefficient: 

g' = a , p^ n " 1 ^ //2 t'" n Sc m (MU)^ n+1 ^' /2 (D2) 

This would involve changes both in a' and n in the case of equation (D2), since, as 
stated in Correlation of Mass Flux With Hydrodynamic Factors, a’ includes numerical 
quantities raised to the power n. One would expect, in the case of equation (Dl), to dis- 
cover that e and n remain constant only if P is properly defined with respect to sur- 
face factors. 
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